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Abstract. We present sensitive interferometric "'^^CO, 
^•^CO and HCN observations of the barred spiral galaxy 
NGC 7479, known to be one of the few barred galaxies 
with a continuous gas-filled bar. We focus on the investi- 
gation and interpretation of -'^^CO/^'^CO line intensity ra- 
tios 7?.i2/i3, which is facilitated by having more than 90% 
of the flux in our interferometer maps. The global (9kpc 
by 2.5 kpc) value of 7^i2/i3 is high at 20 - 40. On smaller 
scales (~ 750 pc), 7?.i2/i3 is found to vary dramatically, 
reaching values > 30 in large parts of the bar, but drop- 
ping to values ~ 5, typical for galactic disks, at a ^'^CO 
condensation in the southern part of the bar. We interpret 
these changes in terms of the relative importance of the 
contribution of a diffuse molecular component, character- 
ized by unbound gas that has a moderate optical depth in 
the "'^^00(1— >0) transition. This component dominates the 
""^^CO along the bar and is also likely to play an important 
role in the center of NGC 7479. In the center, the ^^CO and 
the HCN intensity peaks coincide, while the ^^CO peak is 
slightly offset. This can be explained in terms of high gas 
temperature at the ^^CO peak position. Along the bar, 
the relation between the distribution of ^^CO, ^^CO, dust 
lanes and velocity gradient is complex. A southern ^■^CO 
condensation is found offset from the ^^CO ridge that gen- 
erally coincides with the most prominent dust lanes. It is 
possible that strong ^^CO detections along the bar in- 
dicate quiescent conditions, downstream from the major 
bar shock. Still, these condensations are found close to 
high velocity gradients. In the central region, the veloc- 
ity gradient is traced much more closely by ^^CO than by 

12C0. 
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1. Introduction 

Bars are regarded as an important transport mechanism of 
material toward the central regions of galaxies, fuelling nu- 
clear starbursts and driving galaxy evolution through the 
conce ntrati on of mas s clos e to the nucleus (e.g. Sakamoto 
et al. 199!;, Combes 1990 and references therein). To ac- 
complish this, gas has to flow inward through the bar; 
however, long bars that are filled with dense gas that 
can be traced in CO are rare and thus likely transient. 
Among the few specimens known are NGC 1530 (Rey- 
naud & Downes |1998D , M 100 (Sempere fc G arcia-Burillo 
19971), UGC2855 (Hiittemeister et al. p^99|), NGC 2903 



and NGC 3627 (Regan et al. |1999|). These bars differ sig- 
nificantly in their properties, i.e. velocity field, linewidth 
and derived shock structure. This is somewhat surprising, 
since the gas should in all cases be responding to a strong 
bar potential in a similar way. The dynamics of this re- 
sponse should then have consequences for the properties 
of the gas in the bar and the nucleus, e.g. produce a diffus e 
component unbound from clouds (e.g. Das & Jog 1995 ). 
The evolutionary state of the bar and the degree of central 
mass concentration may be instrumental in regulating the 
condit ions of the gas in the bar (see models by Athanas- 
soula 1992). More specifically, the gas is funnelled from 
outer a;i-orbits to inner a;2-orbits as the bar evolves in 
time (see e.g. the simulations by Friedli & Benz 1993| ). 
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To better understand the gas flow, the star formation 
efficiency along the bar and the feeding of a starburst, it is 
essential to know the physical state of the gas both along 
the bar and in the central structure. An analysis based 
on more molecules than ^^CO, the most abundant species 
after H2, is necessary to investigate gas properties and to 
go beyond the study of the morphological and kinematic 
structure of a system. Achieving this goal requires the 
detection and high-resolution mapping of very faint lines, 
observations which are at the sensitivity limit of today's 
millimeter interferometers. 
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Fig. 1. (a) The map of the ^^CO (1 0) total integrated intensity. The contour levels range from 7.6Kkms~^ 
(1.6 Jy bearn"^ kms~^ or 1.5ct for a velocity width of 400kms~^) to 76Kkms~^ (16 Jybeam"^ kms~^) in steps of 
11.4Kkms~^ (2.4 Jy beam~^ kms~^), 2a) (black), then change to a spacing of 38Kkms~^ (8.0 Jy beam~^ kms~^ , 
8a) up to 380Kkms~^ (80 Jybeam"^ kms~^) (white), (b) Representative spectra from the outer regions of the maps, 
demonstrating that the lines in these regions become narrow, but are clearly detected. The spectra correspond to single 
beams centered on the positions indicated, (c) Dispersion map. The contours range from lOkms^^ to 100 km s^^ in 
steps of lOkms"^. The peak dispersion is llOkms^^. 



In this paper, we present such a study, carried out in 
one of the best known gas-rich barred galaxies. NGC 7479, 
classified a (mild) starburst SBb galaxy containing a 
LINER nucleus, has a gaseous bar of ^ 10 kpc projected 
length (for an assumed distance oi D ~ 32Mpc, see Ta- 
ble|]), which has been repeatedly mapped in ^^CO: Sem- 
pere et al. (1995) present single dish data in the 1^0 
and 2 — > 1 transition. Interferometric 1 — > maps have 
been obtained by Quillen et al. ( |1995| ) and recently by 
Laine et al. (1999). Our results include the most sensitive 
^^CO map to date, but we focus on the relation of the 
"'^^CO emission to the first interferometric ^"^00(1 0) 
and HCN(1 ^ 0) data obtained in NGC 7479, and the 
physical state of the gas we can infer from the comparison 
of these molecules to the ^^CO distribution. Preliminary 



results of this study were given in Aalto et al. (1997a) 



2. Observations 



2.1. Interferometric observations at OVRO 

We have obtained interferometric maps of the 1^0 tran- 
sitions of i^CO, i^CO, HCN and in the inner 1' (9 kpc) 
of NGC 7479 using the Caltech six-element Owens Valley 
Radio Observatory (OVRO) millimeter array. All obser- 
vations were centered on the radio continuum peak. The 
quasar 3C454.5 was observed every ^ 15min as a phase 
and amplitude calibrator. 3C345 was used as an additional 
passband calibrator, and flux calibration was done relative 
to the planets Neptune and Uranus. Calibration was done 
with the standard package developped for OVRO. Contin- 
uum emission with a peak flux of 3 - 4 mJybeam"^ has 
been subtracted from all maps. 

The observations were carried out in the low resolution 
and equatorial configurations of the OVRO array between 
February and June 1996. The naturally weighted synthe- 
sized beam sizes are 4'.'6 x 4'.'25 for ^^CO, 4'.'9 x 4'.'05 for 
"CO and 7'.'2 x 5'.'2 for HCN. At i:>=32Mpc, 1" equals 
155 pc. The positional accuracy of structures in the maps 
is ~ 0'.'5 and thus considerably higher than the resolution. 
Due to the limited uu-coverage, the maps are not sensi- 



Hiittemeister et al.: Gas properties in the bar of NGC 7479 
Table 1. Adopted properties of NGC 7479 



R.A. (1950.0)^> Dec (1950.0)^' Distance Type' 



b) 



Size (blue 



L^\j^ PA disk") PA bar°^ 



23'' 02™ 26.37= 12° 03' 10.6" 32 Mpc SBb W??"^ 



' 45° 



25° 



5° 



a) : Radio continuum peak position (Neff & Hutchings |1992| , VLA A- Array at 1.490 MHz) 

b) : from the UGC catalo gue (Nilson |l973| ) 

c) : from Groosbol (1985) r an gle of inclination 

d) : from Sempere et al. (1995) 

e) : Position angle (north-to-east) of the disk major axis (see e.g. the UGC ca talog ue) and the gaseous bar. The position angle 
of the stellar bar is even closer to that of the disk major axis (Sempere et al. 1995). 



tive to structures larger than ~ 15" for CO and ~ 21" for 
HCN. On-source integration times were 7.2 hr for ^^CO, 
23.5 hr for ^^CO and 11.5 hr for HCN. 

The parameters of the observations are summarized in 
Table |. 

For a beamsize of 4 "4, i.e. ^^CO, a brightness tem- 
perature (Tb) of 1 K corresponds to 0.210 Jy beam~^ at 
a wavelength of 2.6 mm. The spectral resolution of our 
data is 4MHz or 10.04kms^^, and the total velocity range 
covered by the autocorrelator is 1120kms~^, centered on 
i;lsr= 2500 km s^^. For the final presentation, the data 
were smoothed to a velocity resolution of 20 km s^^ (^^CO 
and 13 CO) or 40kms"^ (HCN). 

Imaging was done using the NRAO AIPS package. For 
image display and analysis of the data cubes, the ATNF 
Karma package was also used. 

2.2. Single dish observations at OSO 

Single dish spectra of the 1^0 transitions of ^^CO and 
I'^CO were taken at the Onsala Space Observatory 20 m 
telescope in 1997. The central position of NGC 7479 was 
observed with a beamsize of 33" and system tempera- 
tures of 500K - 600K (12C0) and 200K - 300K (^^CO). 
The spectra were smoothed to a velocity resolution of 
30kms"^ 



3. Results 

3.1. The Distribution of the CO Emission 

3.1.1. Morphology and Central Kinematics 

The ^^CO distribution is displayed in Fig. It extends 
in one continuous structure over the whole length of the 
bar, with a number of secondary peaks along the bar axis. 
The total ^^CO flux in the map is 380 Jykms"^ The dis- 
tribution has a central concentration with a (deconvolved) 
FWHM size of 7"5 x 2"6 (1200 x 400pc at D^i2Mpc). 
This structure holds about 30 % of the total CO flux of 
the bar. The peak position is at a = 23'' 02"" 26.37^ and 
5 = 12° 03' 11 '.'6, coincident to within 1"0 with the ra- 
diocontinuum peak. 

In Fig.|^b, we display some representative spectra from 
the outer parts of our ^^CO map, regions where the confi- 



dence level of the integrated intensity maps (which refers 
to the full velocity width in the center of 400 km s"^) drops 
below 2(T. The lines are clearly detected in these parts 
of the bar, but become very narrow. The same effect is 
evident from the dispersion map (Fig.|l|c): The disper- 
sion (calculated along the line-of-sight, (^Wmom) falls off 
sharply toward the bar edges, but remains at 40kms^^ to 
60kms~i close to the major axis. 

To show that the change in dispersion is not a result 
of projected rotation, we display position- velocity (pv) di- 
agrams along the bar major axis (position angle (PA) 5°, 
Fig.§i) and the disk major axis (PA 25°, Fig.|c). The 
velocity width is not systematically larger along the disk 
major axis, where rotation should have the largest effect. 
Also, the regions of narrow lines are not associated with 
the disk minor axis. The only possible effect of rotation 
is a tentative alignment of the 60kms^^ to 80kms^^ ve- 
locity contours in the ^^CO dispersion map (Fig. Ic and 
Fig. ^) with the disk major axis. 

There is no sign of a 'twin peaks' structure (as coined 
by Kenney et al. 1992| ) perpendicular to the bar, perhaps 
suggesting that there are no X2 antibar orbits, occuring 
within an Inner Lindblad resonance (ILR). However, the 
question of the existence of an ILR in NGC 7479 is con- 



troversial and remains undecided: Quillen et al. (1995) de- 
rived the existence of an ILR close to the nucleus, while in 



the models of Sempere et al. (1995) no ILR exists. Laine et 
al. (1998) argue that the bar perturbation close to the cen- 



ter is so strong that the issue cannot be resolved without 
non-linear orbit analysis. In any case, the ILR would be 
located within 700 pc - 800 pc (4" - 5") from the nucleus; 
thus, we would barely resolve a 'twin peak' structure. 

Another possible signature of an ILR is a ring of ma- 
terial close to its positions, due to the expected pile-up of 
gas in the collision region where xi and X2 orbits inter- 
sect. Such a structure may be indicated in our map. There 
certainly is a clear peak not just in intensity but also in 
velocity dispersion (Fig.|l|c) in the center. This dispersion 
peak fragments into a high velocity system, consisting of a 
number of clumps, as is evident from the pv diagrams we 
show in Fig. ||. This system is confined to within ±3" of the 
center and distinctly different from the material within the 
'main' bar, where the velocity changes much more slowly, 
which is seen in panel (a), a cut along the major bar axis. 
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Table 2. Summary of the interferometric observations of NGC 7479 



Line, v (GHz) 


Date 


Field Center 


Mu-Coverage 


Prim. Beam 


Syn. Beam 


Conf."' 


Tsys 


Noise''' 


'^C0(1 ^ 0) 


2/96 


23'' 02™ 26.4= 


6kA " 46 kA 


63" 


4"6 X 4'.'25 


L,E 


~ 500K 


85 mK 


115.271 


4/96 


12° 03' 10.6" 


(15m - 119m) 




PA -9° 








'=*C0(1 0) 


2/96, 3/96 




5.7kA - 44 kA 


66" 


4 "9 X 4 "05 


L,E 


~ 300K 


23 mK 


110.201 


5/96, 6/96 




(15m - 119m) 




PA -45° 








HCN(1 ^ 0) 


2/96 




4.4kA - 35 kA 


80" 


7"2 X 5'.'2 


L,E 


~ 300K 


12 mK 


88.632 


5/96 




(15m - 119m) 




PA -49° 









o): The array configurations we used at OVRO are labeUed L (low resolution) and E (equatorial) 

&): This is the rms noise per channel (20kms~'^ for CO, 40kms~'^ for HCN). On a mjy beam'^^-scale, the values correspond to 
18 mJybeam-^(i2C0), and 5 mjy beam-i("CO and HCN). 




-18" (south) 0" (north) +18" (SW) -3.6" 0" +3.6" (NE) (SW) -7" 0" +7" (NE) 

Offset along Bar Major Axis (5°) Offset along 50° Axis Offset along Disli Major Axis (20° ) 



Fig. 2. a) Position velocity diagram along the major axis of the ^^CO bar (position angle 5° (north-to-east)). Thus, 
the bar is oriented very close to the north-south direction, b) Position velocity diagram along the axis with position 
angle 50° (i.e. 45° relative to the bar major axis). This is the angle that best shows the central rapidly rotating 
structure, c) Position velocity diagram along the disk major axis, at a position angle of 25° (i.e. 20° relative to the 
bar major axis). For all panels, the contours are multiples of 2a (0.17K or 36 mJybeam^^). 



The central rapidly rotating structure is most easily seen 
at a position angle close to 45° (panel (b)). While there is 
a central peak (at wlsr — 2330 km s^^), this peak does not 
dominate the emission: other features at velocities ranging 
from 2200 km s~^ to 2500 km s~^ are almost as strong and 
spatially separated by 1" — 2". 

A high velocity system like this is a signature feature 
of many barred galaxies. It has been used by Binney et 
al. (1991) to argue the presence of a central bar in our 



Galaxy. The orbits of the gas that can cause this feature 



have been modelled by e.g. Garcfa-Burillo & Guelin (1995) 



for the case of the weakly barred edge-on galaxy NGC 891. 
The st ructu re found in NGC 7479 (already noted by Laine 
et al. 1999 ) also greatly resembles the system seen in 



UGC2855, a strongly barred galaxy where, in contrast 
to NGC 7479, quiescent conditio ns are found for the gas 
along the bar (Hiittemeister et al. [19991) . As in UGC 2855, 
the central high velocity system in NGC 7479 can be in- 
terpreted as a clumpy tilted ring or torus close to the ILR 
(which we would then place at 2" — 3" from the nucleus), 
as a fragmented, tilted, rotating nuclear disk that is freely 
fed from material inflowing through the bar without be- 
ing stopped at an ILR, or even as a dynamically decoupled 
inner bar, as tentatively suggested by Baker ( 2000| ). 



Hiittemeister et al.: Gas properties in the bar of NGC 7479 



5 



12, 



(a) 

"CO on 
'CO 



03 30 



o 
o 

u 

o 



12 03 00 




cut a 
(center) 



cut e 



23 02 27.00 26.00 
RA (1950.0) 



23 02 26.50 26.00 
RA (1950.0) 



Fig. 3. The distribution of the ^^CO emission (contours) superposed on the '^^CO intensity map (greyscale). a) 
Intensity integrated over the full velocity width in the center, i.e. 400 km s""'^. The contours start at 1.4Kkms~^ 
(0.3 Jy beam^^ kms~^ or l.Scr for a width of 80kms~^) and are spaced by 1.85Kkms~"'^. b) The inner region of 
the bar, integrated over 200 km s^^ only to bring out emission from narrow lines. The contours start at 0.9Kkms^^ 
(0.2 Jybeam"^ kms~^, 1.5(7 for a velocity width of 40kms~^) with a spacing of 1.35Kkms"^ The crosses mark the 
positions for which spectra are shown in Fig. 5. Thin gray contours refer to ^^CO. The arrows mark the positions of 
the intensity cuts in Fig.||. 



3.1.2. Comparison with Prior CO Maps 

The general morphological structure and velocity field of 
our map are in excell ent a greement with the results re- 
ported by Laine et al. ( 1999 ). However, the ^^CO emission 
in our map is a continuous bar structure, while the struc- 
ture seen by Laine et al. breaks up into separate clumps 
away from the central condensation, a large part of the 
body of the bar being devoid of emission. As can be ex- 
pected if the two maps are consistent, the position of these 
clumps is in exact agreement with the peaks we see em- 
bedded in more extended emission. 

Comparing the total fluxes contained in the maps, we 
find that our flux is higher by a factor of 1.45 than that 
seen by Laine et al. This discrepancy is caused by our 
map being more sensitive by a factor of 2, due to a 
larger beam size combined with a slightly longer integra- 
tion time. Thus, the larger East- West diameter of the bar 
seen in our map and its continuity are not artifacts of 



smearing by lower resolution, but real effects: Our map 
picks up an additional, more smoothly distributed gas 
component traced by ^^CO. We can compare our total 
flux to the flux seen in the single dish map of Sempere 
et al. ( 1995 ), which covered a larger area than our map. 
However, ^ 80% of the flux is contained the region of the 
OVRO primary beam. We find that we recover more than 
90% of this flux, which is consistent with Laine et al. , who 
report recovering 60% of the single dish flux in their less 
sensitive map. Thus, there is very little, if any, missing 
flux in our CO map and we can analyse the complete 
molecular distribution as traced by ^^CO at high resolu- 
tion. This complete recovery of flux is the main distin- 
guishing feature between our CO map and the interfero- 
metric maps published previously, both by Laine et al. and 
earlier by Quillen et al. (1995). 




Fig. 4. Maps of the onc-dimcnsional velocity dispersion in the inner region where continuous CO emission is detected, 
a) ^■^CO. The contours range from lOkms^^ to lOOkms"^ in steps of lOkms^^. b) ^^CO. The contour levels are at 
lOkms^^, 15kms^^and then range from 20kms^^ to 80kms~^ in steps of lOkms^^. c) Ratio of the ^^CO velocity 
dispersion to the ^^CO velocity dispersion. The contours label ratios of 1.2, 1.4, 1.6, 2, 3 and 5. 
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Fig. 5. Spectra of the central position (taken at the 
position with the highest ^^CO intensity, i.e. a = 
23*^ 02™ 26.37" and S = 12° 03' 11 '.'6 ). The dashed line 
is ^^CO divided by 10 for better comparability, the solid 
line is ^'^CO and the dotted line HCN. All three spectra 
above correspond to a region with the 7 '.'2 x 5 '.'2 size of 
the HCN beam; the ^^CO and ^"^CO data were smoothed 
to the same beamsize as the HCN data. 20 mJybeam""'^ 
correspond to 73 mK (HCN) or 44mK (^^CO). 

3.2. Distribution of^^CO and HCN 

Since ^^CO and HCN usually originate from regions that 
are more compact and confined than ^^CO emitting gas, 
we consider it very likely that we also recover the entire 



flux for those molecules, being limited only by sensitivity, 
not by missing zero-spacings in the interferometer map. 

The distribution of the ^^CO emission is presented in 
Fig. H While Fig. ^a gives the moment map for the full 
central velocity of 400kms~^ (see Fig.||), panel b is in- 
tegrated over only 200kms~^, so that very narrow lines 
are not lost in the noise. Now, a connection between the 
central distribution and the southern feature becomes ap- 
parent. The example spectra shown in Fig.|| demonstrate 
the narrowness of these lines. The total ^^CO flux we de- 
tect is ^ 9 Jykms"^. 

The ^^CO emission is distributed very differently from 
the ^^CO. The central ^^CO peak is significantly offset to 
the north-west from the ^^CO peak (by 2" or 320 pc in 
linear size). 

While the ^^CO distribution in the central concentra- 
tion and south of the center along the bar runs almost 
exactly north-south, the central ^"^CO distribution is best 
described as extending along a position angle of —30°, 
with a curvature toward the south-east (Fig^. The de- 
viation between the ^^CO and ^^CO morphology is even 
clearer in the southern -'^'^CO extension at offsets 7" to 15" 
to the south along the bar. 

No ^'^CO is seen at corresponding offset along in the 
northern part of the the bar. 

The weak feature offset ^ 20" to the north, close to the 
bar end, has no southern counterpart, but is coincident 
with a peak in the ^^CO distribution and thus likely to 
be real. Further peaks at the bar ends are not detected 
with at a sufficient confidence level, since these regions are 
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Fig. 7. The total integrated intensity of HCN, superposed 
on ^^CO (grayscale and thin grey contours). The HCN 
contour levels start at 1 Kkrns"^ (0.42 Jy beani"^ kms""'^ 
or 1.5(t) and increase in 2a steps. The ellipse in the lower 
left corner denotes the HCN beam and demonstrates that 
the HCN emission is not resolved. 



Fig. 6. Offset spectra for ^^CO divided by 10 (dashed) 
and ^'^CO (solid). The positions are marked in Fig^ and 
correspond to regions with the size of the beam of the CO 
observations (~ 5" x 4", see Table 2). For the ^^CO beam 
used here, 20 mJybeam^^ correspond to 97 mK. 

close to the edges of our primary beam, and any emission 
may be affected by the fall-off in sensitivity as well as 
sidelobes of the central peak. One might expect to recover 
""^•^CO again close to the bar edges, due to the crowding 
of elliptical (xi) streamlines, resulting in increased cloud 
collisions and the formation of dense shocked gas. Further 
high sensitivity studies are needed to clarify this point. 

Still, very dramatic changes in the CO to CO line 
intensity ratio are evident, both in the central region and 
along the bar. These changes, which are indicators of 
changing gas properties, will be discussed in detail in the 
subsequent section. 

HCN is clearly detected only in the very center of the 
^^CO distribution (Fig.^. It is unresolved and its position 
is consistent with the ^^CO peak, i.e. offset from the max- 
imum -'^^CO intensity. The total HCN flux in this central 
structure is l.GJykms . 




2000 2500 
Velocity [km/ s] 



Fig. 8. Single dish observations of -^^CO and ^^CO taken 
with the OSO 20m telescope. The beamsize is 33", the 
intensity scale is Tmb (main beam efflciency ?7mb = 0.5). 

3.3. Line intensity ratios 

We present line intensity ratios taken between ^^CO and 
^^CO (7?.i2/i3) in different positions and over different ar- 
eas in Table H In this table, we have computed 7?.i2/i3 and 
its errors individually for the peak intensity, the integrated 
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5 
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Fig. 9. East- West cuts in integrated intensity through 
the bar at the positions marked in Fig.^. The north-south 
width of the cuts correspond to the 4" — 5" beamwidth of 
the CO observations. SoUd Hues are ^^CO, dashed hues 
are ^^CO divided by 20. The integrated intensities re- 
fer to a velocity range of 450kms~^, i.e. the fuU range 
where ^^CO emission is found. The thin dashed hne is 
the 2(7 hmit for "CO and a velocity width (FWHM) of 
40kms~^, i.e. the typical width of off-center ^^CO lines. 
Thus lower ^^CO intensities have no meaning. All offsets 
are relative to the right ascension of the ^^CO peak posi- 
tion. 1 Jybeam"^ kms~^ corresponds to 4.8Kkms~^. 



intensity as measured within the velocity window where 
^^CO is detected and the velocity range within which we 
find ^^CO emission. The respective velocity ranges are 
given along with the values of 7^12/13. Again, it is obvi- 
ous that the ^^CO line usually is narrower than the ^^CO 
line, especially at and close to the ^^CO peak region along 
the bar. Moving away from this peak, the ^^CO linewidth 
drops and approaches that of "'^'^CO (see the dispersion 
maps for ^^CO and ^^CO (Fig.^) and the spectra dis- 
played in Fig. 6). Especially from Fig.^), which displays 
the ratio of the "'^^CO and the ^"^CO velocity dispersion, 
it is clear that the velocity dispersion of ^^CO always is 
smaller than the velocity dispersion of ^'^CO. The disper- 
sions are similar (ratio < 2) in the central 10", where 
the ^'^CO line is fairly wide ((5v,nom « 60 — 80kms^^) 
and in the southern ^^CO concentration, where the ^^CO 
transition becomes narrow (5wmom ~ 10 — 30kms~^). The 
ratio of the dispersion is high (3 - 6) close to the ^^CO 
peak region along the bar, where the '^^CO line is wide, 
while the ^^CO line is narrow. 

The global ^^CO/^^CO intensity ratio, i.e. the ratio we 
find when comparing the total fluxes for the entire map, 
is quite high at ~ 40, certainly much higher than what is 
typical for gas in galactic disks or even 'normal' galactic 
nuclei (see Table |3| for typical values). It is very unlikely 
that this ratio is an artifact of missing flux, since almost 
all flux in ^^CO is seen and we do not expect to miss 
more flux in the the more compactly distributed ^'^CO. 
However, lack of sensitivity to ^^CO is a concern. 

To evaluate the possible magnitude of this effect, we as- 
sume that ^^CO emission in a velocity range of 200 km s^^ 
is present at the 2a level in all places where ^^CO is de- 
tected and we do not find ^^CO. This (unlikely) scenario 
would result in an additional flux of 10.8Jykms~^, i.e. 
roughly double our ^'^CO flux and bring the global ratio 
down from 42 to 19, which we thus regard as a firm lower 
limit. We have also estimated the global value of 7^12/13 
using the global spectra of ^^CO and ^^CO, i.e. spectra 
constructed for the entire spatial extent of ^^CO emission. 
Here, ^^CO is only tentatively detected. Using a 3a limit 
and a velocity range of 450kms~^, we obtain a ratio of 
7^12/13 ^^30, in full agreement with the ratio determined 
using integrated fluxes. 

In the central condensation, the variation in 7?.i2/i3 
goes along with the morphological shift described above. 
7^12/13 ranges from ~ 10 close to the ^'^CO 'ridge' to ~ 30 
at distances of 5" or more from it. The average ratio over 
the central condensation is fully consistent with the ratio 
determined from the single dish spectra taken toward the 
central position (Fig.|8|). Since we expect the central con- 
densation to dominate the single dish flux in the central 
33", this result confirms our prior conclusion that there is 
no missing flux in the interferometry maps. 

The even more drastic variation of 7?.i2/i3 in the bar is 
evident from both Table^ and Fig.^. The latter presents 
another easy way of visualizing the changing line ratios. 
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Table 3. i^co/^^CO 1^0 line ratios (7^l2/l3) for NGC 7479 and typical galactic ratios. The 'cuts' refer to Fig.|[ 
Individual random errors are given in parentheses. The error in the scaling of all ratios due to calibration uncertainties 
is estimated to be ~ 20%. 









Integrated Ratios 




Position 


Peak Ratio 


Region of 


Range"' 


Region of 


Range"' 






"CO line 


[kms-i] 


12 CO line 


[kms-i] 


global (bar and center) 


- 


< 42''' 








^^CO maximum intensity 


27(7) 


19(2) 


320 


25(3) 


440 


^''CO maximum intensity 


16(6) 


18(2) 


320 


23(3) 


440 


cent. cond. (average) 


22(8) 


22(4) 


400 


23(4) 


440 


central region (OSO) 


19(4) 


21(3) 


390 


28(4) 


500 


cent. cond. (individual) 




10 - 


30"' 






Peak of cut b 


27(9) 


25(5) 


220 


29(7) 


400 


^^CO peak of cut c 


19(5) 


17(4) 


60 


40(15) 


340 


OU peak oi cut c (east) 


6.4(1.6) 


10(2) 


105 


14(4) 


200 


12 CO peak of cut d 


11(3) 


10(2) 


80 


28(6) 


200 


"CO peak of cut d 


4.7(2.2) 


5.1(2.4) 


80 


17(3) 


180 


12 CO peak of cut e 


24(10) 


22(5) 


100 


73(37) 


260 


I'^CO peak of cut e 


3.5(1.4) 


3.5(0.6) 


100 


3.5(0.6) 


100 


northern ^''CO peak 


12(4) 


10(2) 


140 


14(3) 


220 




no i^CO detections: 








offset +10" along bar 








>15''' 


220 


offset —17" along bar 








>25''' 


280 


offset —21" along bar 








>25''' 


160 


galactic disk 


~ 6") 








centers of 'normal' galaxies 


13 ±6-''' 









a) The velocity range gives the region where the ^^CO or ^''CO line is detected above the noise. 

b) Under extreme assumptions, this can be brought down to 20 (see text) 

c) Range encountered in the central 10", both in peak and integrated ratios 

d) 3(7 limit c alcula ted over the ra nge w here ^^CO is visible. 

e) Polk et al. [l988|; f) Aaho et al. |l995 



It shows cuts in integrated intensities (over a range of The central gas surface density of NGC 7479 (~ 
450kms~^, chosen to include all -^^CO emission) along 2OOOM0pc~^, using the 'standard' CO to N(Il2) con- 
version factor for comparison purposes only, see discus- 
sion below) is high enough to place the nuclear region in 
the range of IR-bright starburst galaxies (Scoville 1991), 
even though ultraluminous IR galaxies (ULIRGs) can have 
core gas surface densities up to an order of magnitude 
higher (e.g. Bryant & Scoville |1999|, who obtained exceed- 



the bar minor axis at a number of offsets. Here, we focus 
on the southern ^^CO extension. As the ^'^CO emission 
shifts eastward from "'^^CO when we move south along the 
bar, the line ratios change from values around 25 in the 
central condensation to ~ 40 (integrated intensity) at the 
^^CO maximum along the bar. The ratios calculated using 
the peak flux and the velocity range of the ^'^CO emission 
are often smaller than the 7?.i2/i3 referring to the velocity 
range of the ^^CO line since in many places the '^'^CO lines 
are narrower than the ^^CO lines. At the declination of the 
maximum of the southern '^'^CO distribution, 7?.i2/i3 drops 
from 73±37 to 3.5±0.6 (integrated) or 24±10 to 3.5±1.4 
(peak) within 5"! Interestingly, in this position, off the 
^^CO peak, the widths of the two transitions are identical. 
These values of 7^i2/i3 are put into perspective by noting 
that 7?.i2/i3 of 6 - 7 is typical for global ratios in galactic 
spiral arms (Polk et al. 1988), which might already include 
the contribution of some diffuse gas, while Galactic GMCs 
typically show 7^12/13 values of 3-5. An 7?.i2/i3 exceeding 
30 is only seen in a few luminous mergers (Casoli et al. 
199|, Aalto et al. |199lD. 



ing 20000 Mq pc'^ for NGC 2623, Mrk 231 and NGC 6240, 
assuming a standard conversion factor). However, the gas 
surface density in the bar of NGC 7479, where the most 
extreme values of 7?.i2/i3 are reached, is far smaller than 
in luminous starbursts, pointing to a different dominant 
mechanism accounting for the high value of 7?.i2/i3. 

The large changes of 7?.i2/i3 have to indicate a very 
significant change in the properties of the emitting gas, 
scenarios for which will be discussed below 

HCN 1 — > emission is detected only in the nucleus, 
with a 12C0/HCN line ratio « 20. 



4. Discussion 
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4-.1. Gas Masses Based on CO Isotopomers 

The total flux in our map would correspond to an H2 mass 
of 4-10^ Mq «/the 'standard' CO to N(H2) conversion fac- 
tor (SCF = Xco = 2.3 • 1020 cm-2 (Kkms"^)"!, Strong 
et al. 1988) were applicable, which is, however, not the 
case (see discussion below). Instead of giving more masses 
based on the standard conversion factor, which are un- 
likely to be very meaningful, we compare the 'standard' 
mass derived from ^^CO to masses from ^^CO obtained 
with the simple LTE assumptions of optically thin emis- 
sion and a kinetic gas temperature of 20 K. Globally^ i.e. 
avaraged over the entire map, the ^'^CO-derived mass is 
a factor of 16 - 27 (for ^^CO/^^CO abundance ratios be- 
tween 30 and 50) lower than the SCF mass. If the upper 
limit for the ^^CO total flux (see below) is taken, the dis- 
crepancy decreases to 8 - 13, still remaining very substan- 
tial. A mass discrepancy of ~ 10 is also derived for the 
central condensation. In ^^CO peaks along the bar, the 
disagreement between the SCF and the ^'^CO LTE mass 
reaches values exceeding 30, while in the ^^CO maximum 
in the southern bar, it comes down to ^ 2. Of course, these 
varying mass discrepancies simply reflect the changing line 
intensity ratios discussed above. 

The SCF mass stiU exceeds the ^^CO LTE mass by a 
factor 2 - 10 if the kinetic gas temperature is raised from 
20 K to 100 K for the ^^CO - based estimate. Only in the 
^■^CO peak in the bar south of the central condensation the 
SCF mass drops to about half the ^'^CO LTE mass under 
these conditions, which are, of course, not realistic for a 
^■^CO condensation, presumably consisting of dense gas. 
However, it seems that the SCF might be correct in this 
(and only in this) location, which has almost exactly the 
i^CO/i^CO hne ratio found in Galactic GMCs. This can 
be considered as a way to confirm the SCF for a normal 
spiral environment. 

Of course, assuming LTE conditions when deriving 
A^(H2) from ^^CO could underestimate the true_/V(H2) 
when non-LTE conditions apply. Padoan et al. ( 2000 ) find 
a discrepancy of a factor 1.3 - 7 between 'LTE masses' and 
'true' masses in their cloud models. This discrepancy is 
smaller than the difference between the ^'^CO LTE masses 
and the SCF masses we find both for the gas in the center 
and, even more pronounced, for the gas in ^^CO peaks 
along the bar away from center. Still, the ^^CO masses 
derived for LTE should be considered as a lower limit to 
the real mass, while the SCF masses are upper limits. 

Simple non-LTE radiative transfer calculations (using 
the Large Velocity Gradient (LVG) assumption) reproduce 
the results of Padoan et al. to within roughly a factor of 
two and can thus serve as a guide to possible scenarios. 
These calculations show that the ^^00-10-112 conversion 
factor can be brought down significantly (by a factor of 
~ 6) in low density gas (n(Il2) lOOcm"^). Then, high 
i^CO/i^cO 1 line intensity ratios of 20 - 30 are 



predicted. Thus, this scenario points to the presence of 
'diffuse' gas, discussed in detail in the next section. 

When the density is reduced even further, in the LVG 
models the ^2CO-to-H2 conversion factor rises again. The 
SCF can be recovered for a II2 density of ~ lOcm"'^, 
while 7?,i2/i3 remains high. We consider these extremely 
low densities required for the bulk of the gas unlikely for 
the following reasons: To reach a column density that is 
sufficient to ensure shielding of the molecules against UV 
radiation and to be consistent with the observed bright- 
ness of the lines, even for an extreme and very unlikely 
beam filling factor of unity, the models require an unre- 
alistically low velocity dispersion for a given pathlength 
(< 0.05kms^^ pc^^). This is in contradiction to the ex- 
pectation that diffuse gas should be characterized by a 
large velocity dispersion (see below). Another argument 
against extremely low density gas is the expected very 
low i2CO(2-^l)/(1^0) line intensity rati o (< 0.3), in con- 
tradiction to single-dish (Sempere et al. 1995) and recent 
interferometric observations of -'^200(2^1) (Baker 2000). 



4.. 2. Physical Conditions - Diffuse Molecular Gas 

In this section, we present simple scenarios that might 
explain the variations in 7?,i2/i3. The key component of 
these scenarios is the presence of a diffuse or 'intercloud' 
molecular medium (ICM). This medium consists of low 
density (n^lO'^ cm~^), gravitationally unbound, molecu- 
lar gas. For a given column density, iV, the ^■^CO 1 — > 
line emission from gravitationally unbound gas will be 
weaker than that from self-gravitating clouds, because the 
velocity dispersion, 5v, is higher for unbound material and 
thus N/5v is lower for the unbound gas. A low N/5v re- 
sults in a low optical depth and, in low-density gas, also 
results in a low excitation temperature of the transition 
because of reduced radiative trapping. Hence, the ICM 
will be difficult to detect in the ^^CO 1^0 line. One 
may see the ICM as somewhat analogous to high latitude 
molecular clouds detected in the disk of our Galaxy. These 
are also low density structures which, unlike most of the 
disk molecular clouds, are not in virial equilibrium. 

While, under these conditions, ^^CO is optically thin 
and subthermal, ^^CO reaches a moderate optical depth 
close to unity. Thus, "'^^CO still radiates efficiently (and 
is self-shielding) while the ^^CO intensity falls off. This 
moderate optical depth of the ^^CO 1^0 line is the main 
reason why the 'standard' Xqq conversion factor cannot 
be applied to regions where a significant amount of diffuse 
gas is found: the gas is neither optically thick nor virial- 
ized and the 'standard' Xqq overestimates the molecular 
gas mass by up to an order of magnitude. This situation 
is encountered, e.g., in our Galactic center (Dahmen et al. 
1998|). 
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Fig. 10. (a) The integrated ^^CO intensity (white contours) superposed on the velocity field derived from^^CO 
(grayscale and thin gray contours). The velocity contours range from 2280 kms"^ to 2440 kms"^ and increase by 

lOkms"""^. (b) The integrated "'^'^CO intensity (white contours) superposed on the velocity field from panel (a). The 
velocity contours are the same as in panel (a) and thus facilitate a direct comparison of the two panels, (c): ^^CO 
channel map for the velocity range 2280 km s~^ - 2340 km s~^ (white contours) overlaid on the ^^CO velocity field in 

the central part of the bar (grayscale and gray contours). The intensity contours arc multiples of 2(t (5.9Kkms~"'^ 
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4.2.1. The Bar - Outside the Central 10" Diameter 

The large variations in 7?.i2/i3 along the bar, with values 
ranging from 5 to > 40 (and even 73 in one extreme case) 
is best explained by the presence of a large amount of 
diffuse (thin, warm, unbound) ICM gas which dominates 
the ^^CO 1— >0 emission and is very distinctly different 
from the dense molecular clouds preferentially detected 
by ^^CO. In the southern part of the bar, where ^'^CO is 
detected, the difference in linewidth is most dramatic close 
to the ^^CO maxima. Here, the ^^CO lines have a width 
of ~ 200kms"\ compared to - 50kms"^ for ^^CO. This 
explains the significant difference in 7?,i2/i3 depending on 
whether integrated or peak intensities are used (Table I). 
The large linewidths of the ^^CO spectra suggest large ve- 
locity dispersions at the ^^CO maxima along the bar. This 
provides a strong argument in favour of a diffuse medium 
contrary to [^^CO]/[^'^CO] abundance variations. While in 
a 'standard' ^^CO emitting ISM, i.e. a medium where low 
J i^CO hues are optically thick, the [i2C0]/["C0] abun- 
dance ratio is not traced by changes in 7^i2/i3, a change 
in the ^^CO abundance itself could contribute to changes 
in 7^12/13. However, if such a variation reflects changes in 
[-'^^C]/[^^C] with galactocentric radius, it should not ex- 
ceed a factor of two (Wilson & Matteucci 1992, Langer & 
Penzias 1990] ). Even a gradient this small and insufhcicnt 
to explain the observed changes in 7?.i2/i3 is almost impos- 
sible to maint ain in the bar due to efficient mixing (see e.g. 
Friedli et al. 1994 for model calculations, Martin & Roy 



lined above, since the HCN is coincident with the "'^^CO 
maximum and offset from ^'^CO. A situation like this was 
sometimes seen as suggesting abundance changes for the 



1994| and Zaritsky et al. |1994| for empirical studies). Thus 



we consider diffuse gas to provide the only viable expla- 
nation for the large variations in 7?.i2/i3 we find along the 
bar of NGC 7479. 

Of course, ^'^CO may be selectively photodissociated 
in the diffuse medium and thus really be under abundant, 
since it is not self-shielding. However, the efficiency of this 
mechanism also relies on the diffuse nature of the gas. 

According to the single dish map of Sempere et al. 
( |1995D , the 12CO(2-^1)/(1-^0) ratio along the bar is - 
0.6. This is, of course, a global value which is likely to 
vary at higher resolution. It is compatible with cold gas 
at moderate density (10"^ — 10** cm^"^) or warm, thin gas 
(density a few 100 cm~'^). Given the evidence discussed 
aboved, we obviously favour the latter explanation for this 
line ratio. 

The fact that the ^^CO intensity maximum is clearly 
offset along the bar minor axis from the ^■^CO maximum 
in the southern bar indicates that ^'^CO traces a largely 
different component. The relation of the two tracers to 
dust, shocks and star formation will be discussed further 
in Section 4.3. 



4.2.2. The Inner Region - Inside the Central 10" 

At first glance, the situation in the center of NGC 7479 
might be taken to argue against the ICM scenario out- 



CO isotopomers (e.g. Casoli et al. 1992), since HCN re- 
quires high gas densities of at least 10* cm~^ to be sig- 
nificantly excited and may be identified with the region 
where most of the mass resides. In a turbulent environ- 
ment like the nuclear region of NGC 7479, we consider the 
alternative explanation of a kinetic temperature gradient 
between the ^"^ CO and HCN peaks to be more likely. For 
the high gas densities we expect in the central region (i.e. 
n(H2) > lO^cm""^), the higher kinetic temperature would 
occur at the HCN peak rather than the ^'^CO peak be- 
cause the brightness of the ^^CO 1^0 line is a decreasing 
function of kinetic temperature (for temperatures above 
8K in the LTE, optically thin limit). 

The HCN molecule has a high dipole moment and is 
thus much more sensitive to density variations than ^^CO 
and ^■^CO. Thus the HCN peak is not expected to coin- 
cide with the ^^CO peak if the latter is at lower densities 
than the HCN 1— >0 critical density. This is true even if the 
HCN 1^0 transition is optically thin, due to the spreading 
of the molecules over many rotational levels that occurs 
at high densities (i.e. 10^ cm^"^ for HCN) because of the 
high collision rates at these densities. 

We can give an order-of-magnitude estimate of the gas 
density that is likely to occur at the HCN peak: Hav- 
ing the bulk of the gas in the central beam at densities 
of 10^ cm"'^ requires a very low volume filling factor 
of ^ 10^*, which seems unlikely in the center of a galaxy 
undergoing a mild starburst. Also, very high densities sug- 
gest a velocity dispersion per pathlength that is too high 
to be realistic (several lOOkms"^ pc~^). If, on the other 
hand, the density was low enough for HCN to be subther- 
mally excited, the HCN 1— >0 line would be optically thick, 
because the molecules would pile up in the J = level. 
However, in this case (at peak densities of ^ 10"* cm""^) 
we would expect the emission to be more widespread. Low 
densities close to 10'^ cm^'^ lead to very subthermal con- 
ditions and very weak HCN 1— s-0 emission. Thus, we con- 
sider ~ 10^ cm~'^ to be a reasonable density at the HCN 
peak. 

In summary, since the HCN distribution is confined 
closely to the ^^CO peak, we expect this region to also 
be a peak in density (but not necessarily in column den- 
sity), in addition to being a peak in Tkin- In contrast, the 
^^COl^O emission will preferentially trace cool and/or 
lower density gas. 

A kinetic temperature gradient was also used to ex- 
plain the different spatial extents of the ^'^C0(3— >2) and 



HCN emission in NGC 253 (Wall et al. |1991|). The kinetic 
temperature gradient in NGC 7479 leads to the prediction 
that the ^^G0{2^\) / ratio should be higher in the 
HCN peak than in the ^'^CO peak, as is seen in the nucleus 
of IC 694 in the Arp 299 merger (Aalto et al. |l997bD . It also 
means that the nuclear 12CO(2^1)/12cq(^_,o) ratio is 



Hiittemeister et al.: Gas properties in the bar of NGC 7479 



13 



expected to be high, at least close to unity; this seems 
to be the case (Baker 20001 ). ^^CO 2— >1 observations will 
be crucial to further investigate the gas properties in the 
central 1.5 kpc of NGC 7479 and make comparisons to lu- 



minous mergers, where Casoli et al. (1992) and Taniguchi 



et al. (1999), based on data by Aalto et al. |l995| and Casoh 
et al., have argued that ^^CO 2^1 is also relatively faint 
and thus ^'^CO may be depressed, at least in some cases. 

We note, however, ^^CO 1^0 is not as faint at the 
""^^CO maximum as under the extreme conditions of merg- 
ers: Fig.|^ shows the line to be as strong as HCN - it just 
peaks in a different position. The central ^^CO/HCN hue 
ratio of ^ 20 is on the high side if compared to the range 
found in the central region of a number of galaxies, where 
10 is more typical; however the scatter is very large even 



among barred starburst spirals (Contini et al. 1997). Still, 
a value of 20 may be an indication that the amount of 
warm dense gas in the nuclear region of NGC 7479 is only 
moderate and confined to a small region. 

Diffuse gas is likely to be also important even in the 
central condensation of NGC 7479. This is indicated by 
the high values 7^12/13 reaches away from the curved 
'"'^^CO- ridge' in the central condensation. In these areas, 
the "'^^CO lines are much narrower than the ^^CO transi- 
tions (150 km s"^ - 300kms~^ for ^^CO versus 50kms"^ 
- 80kms~^ for ^■^CO). These differences in line shape are 
also obvious from the global single dish spectra (Fig.^). 

Diffuse molecular gas has been detected in the cen- 
ters of a large number of barred and starburst galaxies 
such as IC 342 (Downes et al. p^92| ) or NGC 1808 (Aaho 
et al. 1994) and, of course, in the bar region of the Milky 
Way (Dahmen et al. 1998). It has also been detected in 
the center of the elliptical galaxy NGC 759 (Wiklind et al. 



1997) where a low line ratio of 0.4 is indicative of a two 



component medium made up of dense, cold clumps em- 
bedded in a warm, diffuse molecular medium. A scenario 
where ^^CO and ^'^CO emission arises from separate com- 
ponents was suggested for IC 342 as early as 1990 by Wall 
& Jaffe. The idea that diffuse gas characterized by ^^CO 
emission of moderate optical depth is important even in 
regions of extreme star formation, usually associated with 
large amounts of dense molecular gas, has recently gained 
credibility even for ultraluminous infrared galaxies (Aalto 
et al. [I995i Downes & Solomon |l998D . 

It is interesting to note that the ^'^CO ridge in the 
central condensation is precisely aligned with the steepest 
velocity gradient (derived from ^^CO, Fig. 10). 



4-3. Production of an ICM by tidal disruption and cloud 
collisions 

Cloud disruption due to the tidal field along a bar is one 
mechanism that could produce a diffuse ICM, especially 
close to the central condensation. Another possible source 
of the ICM are off-center cloud collisions. 



GAUXY RADIUS (kpr) 



Fig. 11. Tidal field per unit mass and across unit length 
(see text for explanation) along the bar. The dotted lines 
indicate the tidal field at radii Ikpc and 5.5 kpc, which 
are the semiminor and semimajor axes of the bar. This 
illustrates the change in tidal force that a cloud on an 
elongated bar orbit can experience. 



(i) Cloud evaporation due to the tidal field : As a bound 
cloud moves in a bar potential or elliptical orbit, the tidal 
field across it will vary in time. This produces both inter- 
nal heating of the cloud and clump evaporation from the 
outer regions of the cloud (Das & Jog 1995). The evap- 
orated cloud mass will become part of the low density, 
molecular ICM, raising 7?.i2/i3. 

To get a lower estimate of how the tidal field changes 
over the bar, we first determined the potential of the 
galaxy from the deprojected K image data (Combes, pri- 
vate communication). The potential was azimuthally av- 
eraged to determine ^(r), which was fitted with a poly- 
nomial function that was used to derive a rotation curve 
for the galaxy. The mass to luminosity ratio [M/L) was 
chosen so as to obtain the best fit of the derived rotation 
curve to that observed. We then used this azimuthally av- 
eraged potential to calculate the tidal field over the galaxy. 
Though this does not give the exact variation of the tidal 
field across a cloud in the bar potential, it will give a first 
appoximation as to how the tidal field varies. In Fig[ll| 
we have plotted the tidal field per unit mass, across unit 



length, tf = (dF/dr) (kms~^pc^^)^ against radial dis- 
tance r in kpc. 

With a projected bar length of approximately 10.4 kpc 
and width 1.6 kpc (Fig.|^), an inclination angle of 41° and 
a position angle 25° for the bar (Sempere at al. 1995), the 
deprojected bar length is about 11 kpc and its width is 



about 2 kpc (Arnaboldi et al. 1995). So at the semimajor 
axis distance of r=5.5 kpc, tf = 0.8 ■ 10^"^ (kms~^ pc~^)^ 
and at the semiminor axis distance of r = Ikpc, tf = 



-5.5 • 10 ^ (kms pc 



-1\2 



Thus, the tidal field changes 



from being disruptive at the bar ends to compressive close 
to the bar center. This considerable change in tidal field. 
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Fig. 12. Superposition of the ^^CO distribution (right) 
and the ^"^CO distribution on an archival HST (Planetary 
Camera) red image that shows the shape of the dust lanes. 



which is steepest in the central region, may produce sig- 
nificant cloud evaporation, thus leading to large amounts 
of diffuse gas. 

(ii) Cloud Collisions: Hydrodynamic simulations have 
shown that off-center cloud collisions lead to gas being 
sheared off the colliding clouds, forming trailing extensions 
in the interstellar medium; the gas then dissipates into the 
ISM (Hausman 1981 ). This gas may form part of the dif- 
fuse ICM. In bars, the crowding of closed xi orbits at the 
bar ends enhances the cloud collision frequency signigi- 
cantly and also produces shocked gas. Collisions also cause 
the clouds to lose angular momentum and sink inwards 
leading to the central buildup of gas. This increased con- 
centration of gas in the centers of starburst galaxies and 
barred galaxies means the rate of formation of diffuse gas 
in these galaxies will also be high. NGC 7479 has a strong 
bar and may have undergone a recent minor merger (Laine 
& Heller 1999). Thus, the cloud collision rate should be 
high and there should be a considerable amount of diffuse 
gas in the central ~ 1.5 kpc. 

Distinguishing between the two mechanisms observa- 
tionally is not easy. The diffuse gas produced by cloud 
collisions may often be accompanied by shocked, dense 
gas. Tidally evaporated gas will not have this association. 
However, diffuse gas is expected to spread over the bar 
more quickly than dense gas, which will make it difficult 
to determine its origin. Still, in the future sensitive, high 
resolution observations of a shock tracer like SiO may be 
useful to decide this question. 



4-. 4- The Relation between Tracers of Shocks and Star 
Formation 

The S-shaped, complex velocity field of the gas in the bar 
of NGC 7479 (see Fig p^) has been convincingly explained 
by e.g. Laine et al. ( 1999] ) as the result of gas streaming 
motions in a strong bar potential. They further argue that 
Ha, dust lanes, a large velocity gradient and strong ^^CO 
emission more or less coincide. These are all taken to trace, 
in some sense, shocks, gas compression and star formation. 

Including ^"^CO and the varying line ratios, indicative 
of at least two different components of molecular gas, (and 
seeing all the ^^CO) leads to a more complex picture. In 
Fig. ^, we overlay the ^^CO (right) and ^^CO (left) inten- 
sity distributions on an HST red image retrieved from the 
archive. In the north, where no -'^^CO is detected outside 
the central region, ^^CO follows the strong dust lane al- 
most perfectly. In the south, the visual dust lane is weaker 
or more diffuse and the correlation with the ^^CO distri- 
bution is much less evident. However, the southern ^^CO 
peak coincides with a (faint) peak in a dust lane. South of 
the center, 7?.i2/i3 grows when we move upstream (west) 
across the bar (assuming that the spiral arms are trail- 
ing). In other words, the ^'^CO concentration is located 
downstream from the ^^CO maximum, and possibly also 
downstream from the peak of the shock, if we identify the 
^^CO maximum with the shockfront. This may be justified 
by the close coincidence of the "'^^CO maximum with the 
optical dust lane in the north and also by the excellent 
agreement Laine et al. ( 1999 ) find between ^^CO peaks 
(both in the north and in the south) and the maxima in 
a J — A' NIR color map that is an even better indicator 
of dust than an optical image. 

In the central concentration, the course of the dust 
lanes is not clear. However, here it is the ^'^CQ that very 
clearly follows the steepest velocity gradient (Fig.[l0|b), 
while ^^CO is distributed fairly independently of the ve- 
locity gradient. ^^CO-emitting gas at a velocity between 
2280 km s~^ and 2340 km s^^ traces the velocity gradi- 
ent somewhat more closely than gas at other velocities 
(Fig.[l0|;), but even in this restricted velocity range the 
^^CO emission does not follow the velocity field as clearly 
as the ^'^CO emission does. A number of spectra from the 
central 15" have their emission peak close to 2300 km s^^, 
and a narrow emission component may be associated with 
this peak (faintly visible in Fig. 4). 

The velocity field based on ^'^CO is displayed in 
Fig. ^0|d) . It is even more dramatically S-shaped than the 
^^CO velocity field. We further note that the -^^CO veloc- 
ity field determined by Laine et al. closely resembles our 
^^CO velocity field, and that their ^^CO distribution in 
the central 15" appears more bent than ours. 

This evidence leads us to the following scenario: The 
^^CO emission and a part of the ^^CO emission, especially 
at mid- velocities, traces dense or at least high column den- 
sity gas that follows the velocity field and is compressed 
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in the region with the steepest velocity gradient. We sus- 
pect that this region coincides with the central dust lane, 
which is, however, not clearly visible in the HST image. In 
terms of ^^CO emissivity, this gas component is not dom- 
inant enough to clearly reduce 7<Li2/i3- Also, the dust lane 
may be fairly narrow and unresolved in our data. Thus, 
this component is seen most clearly in the ^^CO data, but 
Laine et al. also trace it in ^^CO since they are not sen- 
sitive to much of the extended emsission. The ^^CO and 
^^CO distributions in our maps deviate because the ^^CO 
emission, especially but not exclusively at non-central ve- 
locities, is due to diffuse gas that is far less confined to 
the shock region in the central 10" — 15" (~ 2kpc). The 
dramatic S-shape in the ^'^CO velocity field (Fig.|l0|d) indi- 
cates that the strong inner bar of NGC 7479 is the feature 
dominating the kinematics of the denser gas in this cen- 
tral region. As discussed earher (Fig. 2), a rotating disk 
or gas on a;2-orbits, perpendicular to the main bar axis, is 
interacting with the bar on the scale of a few hundred pc 
at the very center. 

It is conceivable that the diffuse gas rises to higher al- 
titudes above the plane of the galaxy than the dense gas, 
which may be less turbulent and more strongly dominated 
by the bar potential, which is defined by the stellar poten- 
tial, since most of the dynamical mass is almost certain to 
be in the form of stars. 

Outside the central 2 kpc, the relation between the var- 
ious tracers changes: both the detached northern and the 
southern ^^CO peaks are found close to, but not coincident 
with, regions of steep velocity gradients (Fig.|l^). The dust 
lane (and presumably the region of the strongest shock) 
is now traced well by the ^^CO emission. Possibly, the 
larger overall amount of gas in the central 2 kpc allows a 
diffuse component to spread throughout this region, while 
in the outer bar the gas associated with the shock remains 
diffuse. 

It seems possible that the ^^CO peaks in the outer bar 
indicate conditions that are more quiescent and favourable 
to the formation of bound clouds and ultimately star for- 
mation activity. The low value of 7^i2/i3 in these com- 
plexes suggests that their properties are similar to GMCs 
in the spiral arms of normal galaxies. The cuts presented 
by Sempere et al. (1995) indicate that the Ha emission 
peaks along with the ^'^CO in the southern condensation, 
i.e. is offset from the ^^CO and main dust lane in this area. 
If the Ha traces star formation, this shift away from a re- 
gion that seems to be totally dominated by diffuse gas may 
reflect that star formation occurs wherever the gas density 
is high enough, while the proximity to the shock front may 
play a secondary role. The fact that the ^'^CO condensa- 
tions are close to (offset 1" — 2") strong (^^CO) velocity 
gradients need not be a contradiction to quiescent con- 
ditions. The peaks of the condensations are downstream 
of the strongest gradient by several hundred pc. In addi- 
tion, the interiors of the condensations could be quiescent 
while the condensations themselves are following a strong 



overall large-scale flow. This question, however, requires 
further study on smaller spatial scales. 

Offsets between molecular gas concentrations and dust 
lanes have recently been reported by Rand et al. (1999) 
for a part of a spiral arm in M83. They used ^^CO as 
a tracer, but see only 2% - 5% of the single dish flux in 
their interferometer map. Thus, they probably filter out 
the diffuse component, so their ^^CO data might trace a 
dense, clumpy component similar to what we see in ^'^CO. 

5. Conclusions 

We have mapped the barred spiral galaxy interferometri- 
cally in the J = 1 ^ transitions of ^^CO, ^^CO and 
HCN. Our main results are as follows: 

1. The ^^CO map shows a continuous gas distribution 
all allong the bar. Comparison to single dish obser- 
vation shows that the interferometer does not miss a 
significant amount of flux. The velocity field derived 
from the ^^CO map is complex, showing the S-shaped 
isovelocity contours typical of noncircular gas orbits in 
a strong bar. 

2. A high velocity feature is identified close to the center. 
This may be a ring associated with an Inner Lindblad 
resonance, a tilted rotating disk fed directly by mass 
infall along the bar or even an inner bar decoupled 
from the main bar. 

3. ^'^CO emission is detected in the central condensation, 
the southern part of the bar and a single location in 
the northern part of the bar. HCN emission is only 
detected from the center. The ^'^CO central emission 
is offset by ~ 2" from the ^^CO and HCN intensity 
maxima, which are coincident. Along the bar, the most 
prominent peak of southern -'^'^CO condensation is also 
clearly offset from the ""^^CO distribution. 

4. The ^^CO/^^CO fine intensity ratio, 7?.i2/i3, varies 
dramatically. Globally (9 kpc x 2.5 kpc), 72.i2/i3 is 20 
- 40, a high value compared to typical ratios found 
in the disk component of galaxies or even central re- 
gions of normal galaxies. This indicates a prominent 
contribution of diffuse, unbound molecular gas with a 
moderate optical depth in the ^^C0(1— >0) transition 
in both the bar and the center of NGC 7479. 

5. On smaller sizescales of ~ 750 pc, 7^i2/i3 exceeds 30 
in large parts of the bar, reaching values usually found 
in starburst mergers. Since values as low as 5 are also 
found in the bar, close to the ^'^CO condensation, and 
since a bar environment is very well mixed, we discard 
an underabundance of the ^'^C isotope as a possible 
explanation of the very high 7^12/13 found in many 
places. Instead, this is explained by a dominant com- 
ponent of diffuse gas, readily produced by either tidal 
disruption or cloud collisions in the bar potential. The 
large variation in 7?.i2/i3 is refiected by large changes 
of likely values of the conversion factor from ^^CO in- 
tensity ro H2 column density. In the central 1.5 kpc, 
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the Galactic 'standard' conversion factor (SCF) over- 
estimates the gas mass by a factor of up to 10; in ^^CO 
peaks along the bar the discrepancy is even larger. 
Only in a ^^CO complex in the bar we find the SCF 
to be correct. 

6. The offset in the central HCN (and ^^CO) peak from 

the ^"^CO peak can be attributed to a gradient in ki- 
netic temperature in which the highest gas kinetic tem- 
perature is at the position of the HCN peak. This leads 
to the prediction that the i-''CO(2^1)/(1^0) intensity 
ratio should be higher at the HCN peak than at the 
13 CO peak. 

7. The region along the bar where 7?.i2/i3 is small might 
be an area where the conditions are more quiescent, 
which is also indicated by the narrowness of the lines, 
both in i^CO and ^^CO, found here. If the ^^CO ridge 
along the bar, which coincides closely with the dust 
lanes, is taken to be the location of the bar shock, 
the ^"^CO condensation is behind or downstream this 
shock, possibly in a region where the disrupted (but 
also compressed) gas emerging from the shock can 
form bound molecular complexes. However, the re- 
gion where the velocity gradient is steepest along the 
bar does not coincide exactly with either the ^^CO 
or the ^"^CO distribution. In the center, ^''CO traces 
the steepest velocity gradient much more closely than 
i^CO. Thus, the relation between the molecular trac- 
ers and the shock is complex. 
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